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ABSTRACT 

We present a numerical study of the evolution of galaxy clustering when galaxies flow passively from 
high redshift, respecting the continuity equation throughout. While passive flow is a special case of 
galaxy evolution, it allows a well-defined study of galaxy ancestry and serves as an interesting limit to 
be compared to non-passive cases. We use dissipationless N-body simulations, assign galaxies to massive 
halos aX z ~ 1 and z ~ 2 using various halo occupation distribution (HOD) models, and trace these 
galaxy particles to lower redshift while conserving their number. We find that passive flow results in an 
asymptotic convergence at low redshift in the HOD and in galaxy clustering on scales above ~ 3/i~^ Mpc 
for a wide range of initial HODs. As galaxies become less biased with respect to mass asymptotically 
with time, the HOD parameters evolve such that A/i/Afmin decreases while a converges toward unity, 
where Mmin is the characteristic halo mass to host a central galaxy, Mi is the halo mass to host one 
satellite galaxy, and a is the power-law index in the halo-mass dependence of the average number of 
satellites per halo. The satellite populations converge toward the Poisson distribution at low redshift. 
The convergence is robust for different number densities and is enhanced when galaxies evolve from 
higher redshift. We compare our results with the observed Luminous Red Galaxy (LRG) sample from 
Sloan Digital Sky Survey that has the same number density. We claim that if LRGs have experienced 
a strict passive flow, their (Ag(M)} should be close to a power law with an index of unity in halo mass. 
Discrepancies could be due to dry galaxy merging or new members arising between the initial and the 
final redshifts. The spatial distribution of passively flowing galaxies within halos appears on average 
more concentrated than the halo mass profile at low redshift. The evolution of bias for passively flowing 
galaxies is consistent with linear bias evolution on quasilinear as well as large scales. 

Subject headings: galaxy clustering — passive flow evolution — halo occupation distribution — 
methods: N-body simulations 



1. INTRODUCTION 

An accurate match between galaxies at different red- 
shifts can provide essential clues for constructing theories 
of galaxy evolution. The ancestry of a given population 
of galaxies, however, is complicated to establish: as galax- 
ies undergo major or minor merging with different prob- 
abilities or have different star-formation rates, individual 
galaxies may no longer retain the common internal, obser- 
vational properties at later redshifts. Meanwhile, the evo- 
lution of clustering of the given initial galaxy population 
as a whole can provide another route to trace ancestry, in 
principle, independent of the evolution of the appearance 
(i.e., the internal properties) of galaxies: except for on very 
small scales where the baryonic physics may dominate, the 
clustering of galaxies depends on the clustering of the typ- 
ical underlying density peaks where they reside. Under- 
standing the evolution of clustering of galaxy populations 
can provide a key tool for controlling the uncertainties in 
the evolution of the appearance of galaxies. As a starting 
point; in this paper wc will study how a given galaxy pop- 
ulation is spatially distributed at low redshift under the 
assumption that a single unchanging set of galaxies flows 
from various initial distributions at high redshift. In other 
words, wc consider the effect of the continuity equation in 
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gravitational clustering, with no sources or sink terms due 
to galaxy merging or formation. 

Recent studies show that one can reproduce the clus- 
tering of observed galaxies by taking the halos and their 
subhalos as a proxy between the galaxies and mass with 
a proper correspondence of galaxy luminosity to halo mass 



(IColfn et aLllp gQl'Ncv rinck et al.i26ollKravtsov et al.l2004 
IConrov et all ]2006: Weinberg et all |2006| ). The resulting 
occupancy of the galaxies among halos is often param- 
eterized in the form of halo occupa tion distribution at 
the given redshift (hereafter, HOD) (iBenson et all 120001: 



Peacock & Smith"20U^;'Scliak"2000';'Scoccimarro ct al."2001|; 
Benson 2001: White et al. 2001 : Berlind & Weinberg 2002; 
Berlind et al ll2003l : lKravtsov et a l. 2004; Zch avi et al.ifeoOSbl: 



Zheng et al.l 120051 



Zheng fc Weinberg] i2007i : IZheng et all 



2007al) . In those semi-analytic and SPH simulations, galaxy 



populations are selected by stellar or baryon mass that 
would correspond to observational properties at the given 
time. Such snapshots at different redshifts collectively of- 
fer a general idea on the evolution of clustering and halo 
occupation statistics of galaxies. For exampl e, the distinc- 
tion betwe en old and young p opulations in iBerlind et al.l 
(|2003[ ) and lZheng eFal] (|2005f ) hints an evolution of HOD 
of galaxies with time. 

However, the characteristics of the clustering evolution 
of a given population are difficult to directly infer from 
these studies, as here the set of galaxies at low redshift 
does not necessarily have a one-to-one relation with the set 
at high redshift. Conversely, the connection between pop- 
ulations at different epochs can be tested and confirmed 
once we understand the characteristics of the clustering 
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evolution for a given population by means of N-body sim- 
ulations. 

In addition, once we understand the evolution of cluster- 
ing of galaxies, we also acquire a better handle for the clus- 
tering bias between the galaxies and mass, which is cru- 
cial for deriving cosmological information from the galaxy 
distribution. That is, to isolate useful cosmological infor- 
mation from the observed galaxy distributions, we need a 
good theory to predict the characteristics of galaxy cor- 
relation function for a given population so that we can 
reduce the uncertainty in modeling galaxy clustering. 

What is meant by 'a given population of galaxies' of 
course depends on our definition. By defining a popula- 
tion at the initial redshift regardless of their appearance 
at the observed time, we naturally exclude the addition 
of any new members to the population during evolution. 
The galaxies of a given initial population in reality may 
undergo merging among them, in which case members 
of the population at different redshifts no longer have a 
one-to-one correspondence. In this study, we want to test 
the clustering evolution of a galaxy population when its 
members are preserved throughout evolution. As such, 
the population respects the continuity equation without 
source (i.e., new members) and sink terms (i.e., merging) 
and therefore its number is conserved. We refer to this 
as 'passive flow' evolution of a galaxy population, as an 
analogy to the passive stellar evolution where the stellar 
population within an individual galaxy is preserved. Al- 
though this passive flow is certainly not the whole pic- 
ture of galaxy evolution, it serves as an interesting limit of 
galaxy clustering evolution by itself and so as a template 
to be compared to non-passive flow evolution. It is shown 
that the continuity equation naturally forces the evolution 
of clustering of biased tracers to converge to a small area of 
parameter space, at least on large scales: th e bias on large 
scales converges toward unity with time (|Dekel fc ReesI 



correlation function and HOD that might indicate galax- 
ies entering or leaving the observed population. Hereafter, 
we often use the term galaxy 'evolution' to describe the 
evolution of galaxy distributions through passive flow. We 
will explicitly use 'passive stellar evolution' to distinguish 
from passive flow. An example of such comparison be- 
tween the numerical passive flow evolution and the ob- 
served red galaxies is found in I White et al.l (|2007t ). In this 
paper, we aim at a more extensive study of characteristics 
of passive flow evolution. 

We compare our results with Luminous Red Galaxies 
(hereafter, LRGs) from Sl oan Digital Sky Surve y (here- 
after, SDSS) at z = 0.3 (jEisenstein et al.ll200lh . as an 
observational counterpart. These galaxies reside in very 
massive halos, and their red colors, lack of cold gas, and 
the evolution of luminosity function imply that they have 
had low star formation rates for a few billions of years (e.g., 
IWake et al.|[2006l : lBrown et al.ll2007[ ). The observed small- 
scale clusteri ng implies their low merger rate, at least at 
low redshift (jMasiedi et al.|[2b06l ). Therefore the compari- 
son of our results to the LRGs will manifest signatures of 
non-passive flow evolution in the LRGs if any. With these 
realistic counterparts available, we focus our study of pas- 
sive flow evolution on galaxies that give the LRG number 
density at low redshift. We call these galaxies as 'progen- 
itors' of the LRGs. We do not assume any speciflc obser- 
vational properties of these progenitors, but only assume 
that these progenitors are a population that is distributed 
in halos as a function of the host halo mass at the initial 
redshift and are set to undergo passive flow after the initial 
redshift. We look for any dependence of the outcome on 
the initial redshift. We also study the passive flow evolu- 
tion for lower mass halos to find whether the signatures of 
passive flow depend on a halo mass range. 

In § [2] we describe the parameters for N-body simula- 
tions and methods for assigning LRG progenitors at high 



IQSTtlNusser fc Davisll9"9lllVvlll99dlHamiltonlll998l : lTegmark feeHehtftesI In § [3] we present the results of the first mo- 



19981 ). In this paper we test whether the continuity condi- 
tion, i.e., passive flow, results in further distinct paramet- 
ric signatures. 

We use dissipationless N-body simulations, identify ha- 
los, and assign galaxies to halos based on various halo 
occupation distributions. We trace these galaxies to lower 
redshifts while conserving their numbers, and study the 
evolution of HOD and the correlation function in the turn- 
around regime (i.e., quasilinear or smaller scales). 

We then compare the results of passive flow with ob- 
served galaxy populations. While galaxies flow passively, 
the stellar contents of individual galaxies may or may not 
passively evolve. Even without merging, some galaxies 
within a given initial population may continue to form 
stars with existing cold gas and stay blue, while others 
might end star formation and become red. However, such 
a heterogeneous appearance makes it difficult to select a 
consistent set of objects at the two redshifts, i.e., to avoid 
sources and sink terms. We thus consider red galaxies 
that have lost their cold gas content and ceased star for- 
mation before z ~ 1 or 2 (i.e., galaxies with passive stellar 
evolution) as an observational counterpart. Such galaxies 
will likely maintain a common observational appearance 
during the subsequent evolution. When comparing to the 
observations, we will look for any discrepancies in galaxy 



ments of HOD, the average number of galaxies per halo as 
a function of halo mass, and correlation function of pas- 
sively flowing galaxies. In § 2] we investigate the second 
moments of HOD, the average pair counts within a given 
mass of halo, for passively flowing galaxies. In §[5] we com- 
pare the resulting statistics for passively flowing galaxies 
with the best fit for the observed LRGs, and in § [5] we 
search for a signature of passive flow evolution or environ- 
mental effects in correlation function for the given HOD. 
In §[7] we study the evolution of bias, compared to linear 
theory. In § [5] we discuss the details in the clustering of 
passively flowing galaxies. 

2. SIMULATED HALOS AND GALAXY POPULATIONS 

2.1. N-body simulations 

Our cosmological N-body s imulations use the Hydra code 
(jCouchman. Thomas, fc Pearcc 1995) in coUisionless AP^M 
mode. We use the CMBfast jS eliak & Zaldarriaga 1 9961 : 
IZaldarriaga et al.lll998l : IZaldarriaga fc Seliaki200Q ) linear 
power spectrum to generate many initial Gaussian random 
density flelds at redshift of 49 and evolve them to lower 
redshifts. We generate the initial flelds using the cosmo- 
logical parameters similar to the 1st year Wilkinson Mi- 
crowa ve Anisotropy Probe (WMAP) data (|Spergel et al.l 
[200I : n,n = 0.27, 17a = 0.73, ni, = 0.046, h = 0.72, 
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and n — 0.99. Wc normalize the initial fields by requiring 
(78 = 0-9 at z = and assuming a linear growth func- 
tion. Each simulation represents Vbox = 2h6^h~^ Mpc"^ 
and follows the evolution of 256'^ dark matter particles 
(~ 1.0355 X 10^^ M0 /particle). We compute gravity using 
256'^ force grids with a Plummer softening length of 0.1 
Mpc in comoving unit. A total of 29 simulations are 
used to allow little interference from st atistical variance. 

We use the friends-of-friends method (|Davis "eTallllJSSh 
and identify host halos by adopting a comoving linking 
length of 0.2 Mpc. Wc assign galaxies to dark matter 
halos with various halo occupation models. In our models, 
halos with fewer than 100 particles host at most 0.2% of 
the galaxy population at z = 1 and 3% at z = 2. Wc use 
small group multiplicities only to represent the extreme 
low mass tail of the halo occupation distribution. 

2.2. HOD models and galaxies 

We start with an assertion that the number of galaxies 
in a halo at the initial redshift is only a function of the halo 
mass not of environments. We assume that the initial halo 
occupation distribution follows those observed in the local 
universe and adopt the following form: 

(iVg(M)) = cxp(-A/„in/A/) X [1 + (A//A/i)"] (1) 
where M is the mass of a halo, M^mi is a characteristic 
mass scale for a halo to have one central galaxy, and Mi 
is a mass for a halo with a central galaxy to have one 
satellite galaxy. A central galaxy is assigned to a halo 
based on the nearest integer distribution with the aver- 
age of (iVcon(M)) = exp(-Mmi„/M). Wc find the most 
bound particl(|3 within each halo and label it as a central 
galaxy. For those halos hosting central galaxies, a number 
of satellites are randomly assigned to the rest of particles 
in the halo, based on the Poisson distribution with the 
average of (iVsat(Af)) = (M/Mi)" (that is, the average 
of (7Vsat(M)) = exp(-Afmin/M)(M/Mi)" over all halos). 
The satellite galaxies therefore trace mass inside the halo 
at the initial redshift. 

The observed spatial distribution of galaxies within a 
halo is similar to mass distribution while the subhalo dis- 
tribution from sim ulations is antibiase d with res pect to 
mass and galaxies ( Ghigna et al. 19981. 2000; .Coli'n et al 



1999tlSpringel et al.l200lHDiemand et al.ll2004lGao et al 



2004allbl:lvan den Bosch et al.i2005al:]Nagai fc Kravtsovl2005t 



Zentner et al.l l20d5t [Weinberg et al.l l2006l and observa - 
tional refere nces therein) (but see iTaylor fc Babu]|l2004D 
Meanwhile, iNagai fc Kravtsovl ()2005D and IConrov et al.l 
(|2006( ) showed that tracing subhalos based on their mass at 
the time of accr etion removes most o f the antibias, and us- 
ing this scheme. [Conrov et al.l (|2006t ) reproduced the clus- 
tering of the observed galaxies. Based on these results, 
we assume that satellites of LRG progenitors trace mass 
inside a halo; we cannot be more rigorous, as our mass 
resolution does not allow us to find subhalos or the mass 
of subhalos at any given time. 

We have little information on how the progenitors of 
the LRGs are distributed at z > 1, and therefore we test a 
wide range of initial (A^g(M)). We define 11 different HOD 
models at z = 1 (Models 1-11) and at z = 2 (Models 21- 
31) with variations in Mi/Mmin and a, shown in Table [T] 

*We find a particle with tlie lowest total energy in a halo. 



and [21 at each redshift, 9 models with Afi/Mmin = 2, 10, 
and 25, and a = 0.5, 1, and 2 and additional two models 
with only central galaxies. For Model 10 and Model 30, 
we assign central galaxies without the satellite term in Eq 
[TJ For Model 11 and Model 31, we do not use Eq[T]but 
consider a case in which halos in a very narrow range of 
mass can host central galaxies: we take a mass range of 
a factor 2. We constrain the appropriate initial Mi and 
-A?min by fitting the total number density of galaxies to 
10~^/i^ Mpc^'^, the number density of the observed LRG 
sample with lu minosity limit —23.2 < Mg < —21.2 at 
0.16 < z < 0.36 (jZehavi et al.| [2005al. 

We trace and locate the labeled particles (i.e., galax- 
ies) in the dark matter halos at lower redshifts down to 
z = 0.3 and derive their correlation function. We construct 
HODs of the evolved galaxies by counting the number of 
galaxies per halo at different halo mass bins. We do not 
physically identify which of the galaxies become central 
or satellite galaxies in halos at low redshift, but we make 
a parametric distinction of the central and the satellite 
populations. That is, we estimate {Nccn{M)) by count- 
ing halos with galaxies. When halos host more than one 
galaxy, we count the additional galaxies as satellites and 
calculate {Nsa,t{M)). By construction, the central galax- 
ies at low redshift follow the nearest integer distribution, 
while the probability distribution of satellites will be stud- 
ied in the following sections. From the central and satellite 
number densities we derived, we fit the resulting (iVg(M)) 
to Eq [1] and derive the best fit parameters: Mmin from 
(7Vccn(M)) and Mi and a from (iVsat(M)). The resulting 
(A"g(M)) is not necessarily in the exact form of Eq[l] but 
close enough in most cases that HOD parameters at the 
initial redshift and the final redshift can be compared. Mi 
and a are fitted over halo mass roughly larger than Mi to 
have a better description of the HOD for massive halos. 

2.3. Two-point correlation function 

The mapping between two-point correlation function and 
halo occupation distribution can be described analytically 
by a halo model of two components: a 1-halo term from 
the distribution of excess pairs within the same halo, which 
dominates the small-scale correlation function, and a 2- 
halo term from excess pairs between different halos, which 
dominates the la.rge-scale correlation funct i on (ISeliakll2000l: 



Ma fc Frvll2000l : iBerhnd fc Weinberg|l2002l : ICoorav fc Shethl 
20021) . The distribution of galaxy pairs within and be- 
tween halos depends on the first and the second moments 
of HOD, given the halo mass function and the halo profile. 
Thus features in HOD are closely related to the 1-halo and 
2-halo terms of the correlation function. 

From our simulations, the two-point correlation function 
of galaxies is calculated by counting the number of excess 
pairs at a given separation for each simulation box and 
then averaging spherically over all simulations. As we only 
study the distribution of pairs in this paper, we hereafter 
abbreviate 'two-point correlation function' to 'correlation 
function'. 

3. HOD AND CORRELATION FUNCTION OF PASSIVELY 
FLOWING GALAXIES 

We start by summarizing the characteristics of passive 
flow evolution in the HOD, mainly the first moment (7Vg(M)), 
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Fig. 1. — The evolution of correlation functions and HODs of gala xies that flow passive ly from z = 1. The left panels show ^ in real space 
(upper) and ^ in redshift space (lower). The data points are from [Zehavi et all Il2005al ) for the observed clustering of the LRGs between 
—23.2 < Mg < —21.2 at 2 = 0.3. The middle panels show the evolution of HODs, and the right panels show the decomposed HODs of central 
galaxies and satellite galaxies at 2 = 1 (gray) and z = 0.3 (cyan). Black line: the initial condition at ^ = 1. Blue: at z = 0.3. Red: at 
z = 0.8,0.6,0.4. The input HOD for Model 11, which is a square function, appears smoothed due to our choice of mass bin. The values in 
the parenthesis in the right panels denote {a, Mi/M^i^) at z = 1. 



and the correlation function. The clustering bias of these 
galaxi es will converge toward unity with time on linear 



scales (iDekel fc ReeslllQSTHNusser fc pavisll994l: [Frvll996l: 
iHamiltonl 119981 : iTegmark fc Peebles! 1X9981) . which means 
that the clustering of different models will converge with 
tim^B- In this section, we will find how this translates to 
the way for halos to host galaxies and whether there is 
a similar convergence in correlation function on smaller 
scales. We assume galaxy populations that passively flow 
from z = 1 and z = 2. We use various HOD models that 
span from models with a significant satellite fraction to 
models with little or no satellite galaxies and study the 
difference in their fate. We also find how the result de- 
pends on the number density of galaxies (and hence the 

^The expected bias at z, b^, for passively flowing galaxies from 
zo, is fez = {bzQ — 1)/G-|-1 on linear scales where b^g is the galaxy 
bias at zg and G is the growth factor between z and zq . See § [7] 



mass scale of the halos) . 

3.1. Galaxies flowing passively from z = 1 

We assign galaxies to dark matter halos at z — 1 us- 
ing 11 different HOD models (Models 1-11) and trace the 
resulting evolution of the HODs and correlation functions 
down to redshifts of 0.8, 0.6, 0.4, and 0.3. 

Table [1] lists the HOD parameters at z = 1 and at 
z ~ 0.3 for Models 1-11. Figure [T] shows the evolution 
of correlation functions and HODs for Models 1, 3, 5, 7, 
10, 11 among those listed in Table [1] The black solid 
lines represent the initial HODs and correlation functions 
at z ~ 1, and the blue solid lines are for the resulting 
HODs and correlation functions at z = 0.3 while the red- 
dish solid lines show the values at intermedi ate redshifts. 
The da ta points shown here are taken from IZehavi et al.l 
(|2005aj) for the observed clustering of the LRGs between 
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Fig. 1. — continued 



-23.2 < Mg < -21.2 at z = 0.3, which sets the fiducial 
number density of LRG progenitors in this paper. As ex- 
pected, the galaxy clustering (except for Model 3) grows 
slowly with time on all scales while the bias with respect 
to the dark matter decreases. 

In Models 1, 2, and 3, in which Mi/Mj^in ~ 2, the frac- 
tion of satellite galaxies is much larger compared to those 
with larger values of Mi/M,nin. As massive halos host 
more satellites, low-mass halos host a smaller number of 
central galaxies at a fixed number density. Thus the clus- 
tering is more weighted by high mass halos. The 1-halo 
term is prominent in these models because the large satel- 
lite fraction of massive halos emphasizes the structure of 
a halo. The trend becomes extreme in Model 3 in which 
massive halos are weighted most heavily. The evolution of 
clustering in Model 3 is not as evident as in other models, 
first because of the statistical noise due to the small num- 
ber of very large halos and and second because of the large 
initial bias which drives a decrease in bias toward unity 
faster: as the galaxy clustering converges toward that of 
dark matter, the growth of galaxy clustering with a larger 



initial bias slows down more relative to the growth of dark 
matter clustering. As halos accrete mass and merge, the fi- 
nal (A^g(M)) at z = 0.3 is parameterized with bigger values 
of Mmin and Mi than at z = 1 (Table [1]). The evolution 
of a shows an increase in Model 1, but changes very little 
in Models 2 and 3. In all cases, the fraction of satellites 
increases with time due to halo merging events. 

Models 4, 5, and 6 start with Mi/Afmin 10 at z = 
1 with larger values of Mi and smaller values of Mmin, 
compared to Models 1-3. As a result, the clustering both 
from 1-halo and 2-halo terms is lower than in Models 1-3. 
While Minin increases with time, Mi changes little. The 
value of a increases or decreases toward near unity. Due to 
the steady Mi, {Ng{M)) along different redshifts overlaps 
near Mi. Model 5, with an initial a = 1, shows little 
evolution in a. 

In Models 7, 8, and 9, Mi/M^in - 30 initially. The 
clustering strength is even lower in these models as the 
values of Mn^n is even smaller. As the satellite population 
decreases from Models 7 to 9, the 1-halo term in the corre- 
lation functions is suppressed due to the exclusion effect in 
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Fig. 2.— The initial HODs at 2 = 1 (left) and the final HODs at z = 0.3 (right) for Models 1-11. Black : Models 1-3 (from Solid, 
long-dashed to short-dashed in order). Red : Models 4-6. Blue : Models 7-9. Cyan : Model 10. Green : Model 11. One finds that different 
initial (Afg(M))s asymptotically converge at 2 = 0.3. 




Fig. 3. — The initial ^ at 2 = 1 (left) and the final § at 2 = 0.3 (right) for Models 1-11. Top : ^ in real space. Bottom : § in redshift space. 
Black : Models 1—3 (from Solid, long-dashed to short-dashed in order). Red : Models 4—6. Blue : Models 7-9. Cyan : Model 10. Green : 
Model 11. 
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a halo finder. However, the evolution of clustering quickly 
recovers the exclusion effect at later redshifts. The value 
of Ml decreases with time while Mmin increases with time. 
Again, a evolves toward near unity. 

In Model 10, we initially populate only the halos with 
central galaxies without satellites. At low redshift, many 
of these central galaxies turn into satellite galaxies by halo 
merging, and the results at z = 0.3 are very similar to 
those of Model 9. 

The evolution of HOD for Models 1-10 is therefore char- 
acterized as following: a approaches toward an attractor 
near unity, and a large initial A/i/Af„iin decreases with 
time (due to halo merging without galaxy merging) while 
a small initial Mi/Mmin stalls. Such evolution of HOD 
parameters drives the convergence of (A'g(M)) at low red- 
shift. Figure O and [3] show the comparisons of HOD and 
correlation functions between different models at the ini- 
tial redshift {z = 1) and the final redshift (z = 0.3). 
From the figures, it is clear that passive flow evolution 
leads to an asymptotic convergence of clustering on scales 
above r ~ 3/i~^ Mpc and of the first moment of HOD. 
Most of the models produce a stronger clustering than 
that of the observed LRG data while Model 10 marginally 
fits the data. From Figure [31 clustering on small scales 
(r < 3ft.~^ Mpc) varies considerably despite the similarity 
in their HODs. The difference on small scales implies that 
the second moment of HOD, that is, the average number of 
galaxy pairs within a halo, may differ for different models. 
The difference in redshift-space clustering at z = 0.3 ex- 
tends to the larger scale than in real-space clustering due 
to the finger-of-God effect. We revisit this issue of second 
moments in § |4l 

In all models, the increase in il/,„i,i is consistent with 
the typical mass accretion between z = 1 and z = 0.3 fo r 
the corresponding mass range re.g.. lWechsler et al.|[20M) . 
On the other hand, the mass accretion of halos near Mi 
will often involve accretions of new galaxies into the halos. 
For Models 1-3, the galaxy accretion near Mi is minor 
compared to the mass accretion, as the initial value of Mi 
at z = 1 is closer to Mmin, driving the increase in Mi with 
time. For Models 7-9, the large Mi/Mmin at z = 1 implies 
that the mass accretion to a halo near Mi will often ac- 
company galaxy accretion, which must be efficient enough 
to decrease Mi with time. Models 4-6 shows an interme- 
diate behavior. 

The small values of Mi/Afmin at low redshift, that is, the 
lack of shoulder between Mi and Mmin, is different from 
the predictions of {N^ jM)) for subhalos or galaxies in non- 
passive studies (e.g.. iBerlind et aL 2003t iKravtsov et al.l 
[200llConrov et al.ll2006l: [Weinberg et al.li2006n i. The dis- 
crepancy is a natural result for passive fiow evolution as 
we do not include any merging or tidal disruption of trac- 
ers. With tidal stripping and merging between subhalos 
or merging to the center, the number of subhalos will de- 
crease with time for a broad range of host halo mass. Semi- 
analytic studies and N-body simulations show that the 
resulting cumulative number density of subhalos reaches 
unity at a small subha lo-to-host halo mass ratio (or circu 



This small ratio should roughly represent a large value of 
Mi/Mtnin by the definition of Mmin and Mi. Conversely, 
the cumulative subhalo mas s function wit h no tida l disrup- 



tion an d merging (e.g., Z cntner fc Bullock>2 003: van den Bosch et al 
l2005b[ ) hints small values of Mi/Mmin for our passively 
flowing case. Thus, one can reason that, with non-passive 
flow. Ml would increase, and then, for fixed number den- 
sity, Mmin would have to decrease to include new central 
galaxies (i.e., host halos that have accreted enough mass 
to pass the threshold). With decreasing Mmin (i.e., de- 
creasing mass threshold for subhalos). Mi will be read- 
justed. However the value of Mi/Afmin will remain large, 
as the subhalo mass function is roughly self-similar for 
different host halo mass, although sli ghtly more abun- 



dant in more massive host halos 



e^ 



Gao et al.l[2004bl: 



Ivan den Bosch et al.ll2005bl : IZentner et all 120051 ) 

For galaxies, tidal disruption will be much less efficient 
than for subhalos as they are more tightly bound systems. 
Nevertheless, dynamical friction will eventually merge some 
of them to the central object, more effectively for satel- 
lites s itting in more massive subhalos (|Binnev fc Tremaind 
Il987f ). which will result in a larger Mi/Mmin than our pas- 
sively flowing case. The merging products corresponding 
to the discrepancy between passively and non-passively 
flowing galaxies should have consequences elsewhere, such 
as in LRG luminosity function, so that we can verify the 
model. 

As a caveat, we note that our passively flowing sam- 
ples at the observed time might not simply correspond to 
galaxy populations in non-passive studies, which are con- 
structed by relating the luminosity of galaxies to the sub- 
halos and host halos above a given mass (or velocity) cut, 
whether the selection of the subhalos is based on values at 
the time of accretion or at the observed time. Again, our 
samples are defined at the initial time, not at the observed 
time. 

All the models we discussed so far show the composite 
information on galaxy statistics from halos of various mass. 
To probe the response of the delta function in host-halo 
mass to passive flow evolution, we select a mass range with 
a width 2 which gives the fiducial number density and 
assign one central galaxy per each halo (Model 11). This 
scheme populates more galaxies in lower mass halos than 
the previous models do, depriving massive halos of galaxies 
at the initial redshift. Due to the deprivation of galaxy- 
hosting halos of mass between Mmin and Mi, relative to 
the other models, this model mimics the observed shoulder 
in HOD of the LRGs at z = 0.3. The resulting correlation 
function at z — 0.3 appears consistent with the observed 
data for LRGs as well. The response of a delta function at 
different mass will slightly vary; otherwise, a superposition 
of delta functions of different halo-mass, when weighted 
by the initial {Ng{M)), will correspond to the evolution of 
galaxies in Model 10. The correlation function of Model 
11 does not appear different from other models as much 
as its HOD does. 

3.2. Galaxies flowing passively from z = 2 
We next observe the response of correlation function to 



^e l ocitvratio) ^ K^gentner fc Bullock | 200| Taylor fc Babul^ ^^^.^^ ^^.^^ ^.^^^ ^^^^^^^ ^ ^^^^^ 1^^.^^ 

200i Oguri fc Lei [200^: |Di.emand et al.||200^: |Gao.e^ . = 9 w MnHok 9.t^^ and mn^nam^ho ro.uK^u. o^mh.- 



2004bl : Ivan den Bosch et al.ri2005bt IZentner et all |2005[) 



z = 2 for Models 21-31 and compare the resulting evolu- 
tion of correlation function and HODs at z ~ 0.3 (Figure 
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parenthesis in the right panels denote {a, Mi/Mmin) at z = 2. 



Hand Table HI) with those from z = 1 (Models 1-11). We 
investigate whether and where the convergence of HOD 
parameters occurs if galaxies flow passively from z = 2, 
relative to the galaxies from z = 1. We will show that the 
evolution from z = 2 to z = 0.3 is qualitatively similar 
to the evolution from z = 1 to z = 0.3 except that the 
consecutive time steps cover a broader range of evolution. 

At higher redshift, the typical halos have a lower mass. 
The average mass of host halos for the fiducial number 
density is smaller by more than a factor of 2 than that 
at z = 1, but these halos are more biased than the typ- 
ical host halos at z = 1 due to the rari t y of the halos 
at hi gh redshift ([Press fc Schechterl Il974l : iBardeen et all 
Il986[) . The initial clustering strengths at z = 2 are sim- 
ilar to those from z = 1. and so the resulting correlation 
function at z = 0.3 is larger for Models 21-31 than Models 
1-11. 



In general, the transition from the 1-halo to the 2-halo 
term in the correlation function occurs at a smaller sepa- 
ration at the initial redshift, as expected from the smaller 
halo-size of the dominant halo population and the steeper 
halo mass function at z = 2. Models 21-31 start with a 
similar or slightly smaller initial ratio of satellite to cen- 
tral galaxies compared to Models 1-11 except for Model 
23 and Model 3. However, Models 21-31 produce a larger 
satellite fraction at z = 0.3, which is due to the longer 
time available for halo merging and accretion events. 

In Model 21, 22, and 23, in which Afi/Mmin ~ 2, the 
fraction of satellite galaxies is much larger compared to 
the other models, like in Models 1-3. Both Mmin and Afi 
increase at z = 0.3, and a approaches somewhere near 
unity. 

Model 24, 25, and 26 have initial Mi/A/min ~ 10. Both 
Mmin and Mi increase with time, unlike Models 4-6, but 
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Fig. 4. — continued 



the ratio of Mi/Mmin decreases considerably at z = 0.3, 
as in Models 4-6. Again, the value of a converges toward 
near unity. 

In Model 27, 28, and 29 with Mi/M^i^ ~ 25 initially, 
we observe the strong exclusion effect at z = 2 that quickly 
recovers with time. In these models. Mi at z = 0.3 is less 
than All at z = 2 just like in the corresponding models at 
z = 1 (Models 7-9). Again, a approaches toward unity. 

Model 30 is almost identical to Model 29. Model 31 
evolves to have a region of shoulder at z = 1 but then 
takes on a power-law shape at z ~ 0.3 and produces a 
stronger clustering than the observed LRGs. 

From the evolution of the real-space correlation, wc find 
a feature prevailing in most of the models: the growth 
of clustering is impeded near r ~ 2 - 5/i-i Mpc (but 
r > 3/1-1 jyjpj, fQj. Model 23). That is, in Figured the 
correlation functions at different redshift become squeezed 
together over r ^ 2 — 5h^^ Mpc, implying a suppression 
in growth of galaxy clustering. This could be a signature 
associated with a turn-around and infall in the structure 
formation. In redshift space, the features arc less obvious 



but still traceable. For Models 1-11, the corresponding 
feature is observable on slightly larger scales. We further 
discuss about this feature in § [5] 

Figure O and [5] show the comparisons of HOD and cor- 
relation function for Models 21-31 at the initial redshift 
(z = 2) and the final redshift (z = 0.3). Figure [5] shows a 
stronger convergence of these final HODs than observed 
in Models 1-11. The asymptotic a at z = 0.3 moves 
toward slightly larger a than in Models 1-11 (Figure [51). 
Likewise, Figure [H] shows a stronger convergence between 
models in the large-scale clustering (r > 3h~^ Mpc) at 
z = 0.3 as well as in small-scale clustering, which is due to 
the larger bias at the initial redshift and longer evolution 
time. Again, the small-scale correlation function remains 
different between models although their HODs look nearly 
identical. 

Figure [7] provides a clear view of similarities and dif- 
ferences in the evolution of HOD parameters for the two 
formation redshifts. The middle panel shows that Mmin 
increases with time while a converges toward near unity. 
The range of resulting A/min at z = 0.3 for Models 21- 
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Fig. 6. — The initial ^ at 2 = 2 (left) and the flnal 5 at 2 = 0.3 (right) for Models 21-31. Top : 5 in real space. Bottom : ^ in redshift 
space. Black : Models 21-23 (from Solid, long-dashed to short-dashed in order). Red : Models 24-26. Blue : Models 27-29. Cyan ; Model 
30. Green : Model 31. The convergence in correlation function is stronger when galaxies flow passively from 2 = 2 than from 2 = 1. 
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31 is similar to the range of Afmin at z = 0.3 for Models 
1-11. The right panel shows that Mi tends to converge 
to Ml 2 - 5 X IO^Mq at low redshift while the con- 
vergence is stronger and at lower Mi for Models 21-31. 
The left panel shows that the resulting Mi/Mmin con- 
verges toward Mi/Mmin ~ 3 — 4 while the convergence 
is stronger and at lower value of Mi/Afmin for Models 21- 
31. It appears that the convergence of a is stronger and at 
a slightly larger value, somewhere above unity, for Models 
21-31 than Models 1-11. 

3.3. A test with host halos of lower mass: increasing the 
number density 

We have shown the characteristics of passively flowing 
galaxies for high mass halos that are consistent with the 
number density of the LRGs between —23.2 < Mg < 
—21.2 at z = 0.3. The characteristics could depend on 
the peculiarity (or rarity) of an extreme tail of nonlinear- 
ity, and thus we investigate the effect of passive flow for 
the galaxies of a larger number density, i.e., in lower mass 
halos. As both Models 1-11 and Models 21-31 produced 
a clustering stronger than one observed for the LRGs, we 
also ask whether the observed LRG clustering can be bet- 
ter explained as a random fraction of a parent population 
with a larger number density that has evolved passively to 
lower redshift. 

We assign galaxies to dark matter halos at z = 1 using 
values of Mi/Mmin and a similar to Models 7-10 but with 
/owrtimes the fiducial number density 4x 10~^/i'^ Mpc~^) 
and trace the resulting evolution of (A'g(M)) and two-point 
correlation functions down to redshifts of 0.8, 0.6, 0.4, and 
0.3 (TableEl). We label these models as Models 7n4-10n4. 
These models have slightly larger initial and final satellite 
fractions compared to Models 7-10. 

The evolved clustering of Models 7n4-10n4 is weaker 
at z = 0.3 than the clustering of the LRGs, which im- 
plies that the parent population of the observed LRGs 
will have a number density between 10~'^h^ Mpc""^ and 
4 X 10~'^h^ Mpc""^ under passive flow evolution. The evo- 
lution of Mi/Mniin and a for these models is very similar 
to that for Models 1-11, although the convergence of Mmin 
and Ml is at lower mass than in Models 1-11. Therefore 
the characteristics of passive flow evolution in HOD pa- 



rameters we find, that is, the decrease in A/i/A/min and a, 
are fairly robust for a wide range of halo mass. 

The convergence of Mi /Mj^in and a can be understood 
from the idea that the passive flow drives the galaxies 
to flow and distribute like mass, as gravity does not dis- 
tinguish galaxies from other mass components. The dis- 
tribution of passively flowing galaxies will then converge 
to (iV(M)) of mass, which is a power law with an index 
of unity. With our parameterization, the convergence of 
(7Vg(M)) to the linear equation in halo mass forces a de- 
crease in Mi/Mmin and convergence of a toward slightly 
larger than unity due to our division to central and satellite 
populations. This asymptotic convergence will be less effi- 
cient in mass components lower than Afmin, as host halos 
build up hierarchically to larger ones with time rather than 
disassemble in field, therefore maintaining the biased clus- 
tering. The normalization of (A'g(M)) will be determined 
by distributing the total number of galaxies to halos of 

M > M„in. 

Note that HODs of all models except for Model 11 re- 
semble the old galaxy population defined at the observed 
redshifts (i.e., without passive-flow restricti on) in studies 
of SPH simulations or semi-analytic models (jBerlind et al.l 
[2003,; Zheng et al. 200a) . Model 11 on the other hand re- 
sembles the young galaxy populations in those studies. 

4. THE EVOLUTION OF PAIR COUNTS WITHIN A HALO 

In the previous sections, wc found that while various ini- 
tial models evolve to a well-defined region of (A'g(M)) at 
z = 0.3, there still remains a considerable difference in the 
small-scale clustering at z = 0.3 between models. While 
the clustering on large scales can be modeled by the 2-halo 
term that depends on the excess number of pairs between 
halos, i.e., {Ng{M)) weighted by halo bias and halo mass 
function, the small-scale clustering (1-halo term) depends 
on the average counts of excess pairs within a halo, i.e., 
(Ng{Ng — 1))m, the spatial distribution of the pairs within 
the halo, and halo mass function. Therefore the compar- 
ison of pair counts within a halo will help to explain the 
difference in the small-scale correlation functions. 

When samples are defined at the observed redshifts, 
SPH and semi-analytics calculations imply that the satel- 
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lite probability d istribution is mode led well b y the Pois- 
son distribution JWhite et al.ll2001t rBe rlind fc Weinberg! 
I2002t iKravtsov et al.ll2004 IZheng et al.lf2005[) that gives 
{Ng{Ng - 1))m = (A^g(M))2, while the central galaxies 
follow {Ng{Ng — 1))m = 0. From the pair counts of the 
passively flowing galaxies, we thus will find the effect of 
removing processes, such as merging, destruction, or cre- 
ation of galaxies, on the second moments of HOD, again, 
provided that the initial populations follow HOD as a func- 
tion of halo mass only and have Poisson-distributed initial 
satellites if any. 

Figure [8] shows the pair counts aX z ^ 0.3 for galaxies 
that passively flowing from z = 2 (Models 21, 23, 25, 27, 
30, and 31) and the corresponding small-scale correlation 
function in real space. In the second field of the pan- 
els, the pair counts are divided by (7Vg(M))^ to be com- 
pared to the nearest integer and the Poisson distribution: 
{Ng(Ng - l))/(iVg)2 = for the nearest integer distri- 



bution and unity for the Poisson distribution. For HOD 
models with a negligible satellite fraction for halos at the 
low-mass end where (7Vce„) < 1, (A^g(^g - l))/{Ngf will 
start from zero and approach unity as the satellite fraction 
grows with M (such as Models 26, 28, and 29). On the 
other hand, some of our models allow a non-zero satellite 
fraction for (iVcon) < 1 at the initial redshifts, and this 
reverses the shape of {Ng{Ng — 1))/ (A^g)^ at the low-mass 
end as shown in Figure [5] for Models 1 and 7. 

From the figure, the evolved {Ng{Ng - l))/{Ng)'^ at 
z = 0.3 nicely converges to unity for massive halos in 
all models. This implies that the passive flow evolution 
results in satellite populations that follow a Poisson distri- 
bution over a wide range of halo mass, although some of 
the models, such as Models 29, 30, and 31, started with a 
nearly or purely nearest integer distribution at ^ = 2. In 
low-mass halos, galaxies at z = 0.3 follow a nearest integer 
distribution by construction (§ 12. 2p as long as each halo 
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hosts one or zero central galaxy. The third field of the 
panels, where we show the pair counts of satellites alone 
calculated among halos with central galaxies, reveals fur- 
ther details on satellite distributions. Except for Model 31, 
satellites in models with a smaller initial satellite fraction 
or a larger increase in the satellite fraction between z = 2 
and z = 0.3 (i.e., Model 26, 28, 29, 30), converge better 
to the Poisson distribution, though slightly sub-Poisson. 
Satellites in models with a larger initial satellite fraction 
and those in Model 31 tend to be slightly super-Poisson, 
while converging toward a Poisson distribution at the mas- 
sive end. 

Considering the Poisson distribution of subhalos or satel- 
lite galaxies from non-passive evolution (IWhite et al.ll200ll; 



Berlind fc WeinberEj[200^ : lKravtsov et al.ll2004l : IZheng et all 
20051) . our results implies that the processes of merging, 
formation, or destruction of galaxies (or subhalos), which 
are not included in our passive flow evolution schemes, pre- 
serve the Poisson distribution. This implies that the con- 
vergence toward the Poisson distribution probably arises 
from the random nature of the halo merging (Zheng Zheng, 
a private communication). Although not shown in the 
figure. Models 1-10 show similar results while the conver- 
gence toward the Poisson statistics is not as strong as Mod- 
els 21-31. In Model 11, all galaxies for M > lO^'^M© ap- 
proach a Poisson distribution (i.e., {N^{Ng — 1))/{N^)^ ^ 
1), while being slightly super-Poisson, due to the relatively 
large satellite-to-central galaxy ratio for intermediate mass 
halos despite the small total satellite fraction. Model 31 
has a more moderate transition from central to satellite 
domination. 

The fourth and the fifth fields in Figure [5] show ratios of 
{Ng{Ng - 1))m and (NgCM))"^ with respect to the values 
for Model 30. When we compare the first and the fourth 
fields of the figure, we find that the trend of {Ng{Ng — 1))m 
for M > IO^^Mq is qualitatively consistent with the trend 
of ^ on small scales, even without considering the details 
of the mass-dependent halo profile. The comparison be- 



tween {Ng{Ng - 1))a/ and (iVg(M))2 shows that the differ- 
ence in {Ng(Ng — 1))m is rooted in the small difference in 
(7Vg(M))^ between models (see § (5]). The overall conver- 
gence of {Ng{M)) gives the convergence in clustering, and 
it is stronger on large scales because the 2-halo term de- 
pends on the integrated effect of (7Vg(M)) over halo mass. 
The small-scale clustering is more sensitive to the details of 
(iVg(M)) through {Ng{Ng — 1))m, as each mass bin corre- 
sponds to a range of scale r in ^(r) consistent with a typical 
halo size of that mass. We conclude that the difference in 
small-scale clustering between models is due to the differ- 
ence in (A^g(M)) (and its effect on (7Vg(Afg — 1))m), and not 
to the nature of the scatter (e.g., Poisson or non-Poisson). 

5. COMPARISONS TO THE CURRENT OBSERVATION OF 
THE LRGS 

Previous sections show that passive flow leads to a small 
value of Mi/A/niin and a: in general, Mi/Mmin < 10 and 
a slightly above unity. In this section, we compare our 
HOD parameters for the passive flow with the fit to the 
observed LRGs bv lZheng et all l|2007bl ). Even for the same 
cosmology, the details of the group finding method or the 
resolution of simulation may alter mass scale or mass func- 
tion of halos. Therefore, in order to discuss th e differe nces 
between our results and the fit bv lZheng et al.l ()2007bD . we 
first need to identify an LRG population at z = 0.3 in our 
own simulations that is consistent with the ob served clus- 
tering a nd the shape of (7Vg(M)) adopted bv lZheng et al.l 
(|2007bD . By this, we can calibrate our mass scale and as- 
sess the physical significance of the parameter space of the 
HOD confined by passive flow evolution. 

A five-parameter HOD fit (eq. [2]) to the observed LRG 
clustering, derived for our fiducial cosmology, is kindly pro- 
vided by Zheng Zheng (in the footnote of Tabled]). The 
five-parameter HOD is defined as 

(iVg(M)) = 0.5 [1 + erf[logio(M/M:„iJ/cTM]] 

X [1 + [(M - Mo)/Af{]"] (2) 
where M^^^ is the characteristic minimum mass to host 
a central galaxy, M[ is a mass for a halo with a central 
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galaxy to host one satellite when Mq ^ M[, Mq is the 
truncation mass for satellites, and ct m is the characteristic 
transition width (jZheng et al.ll2005[ ). 

Co mpared to the Jenkins mass function (jJenkins et al.l 
[200l[) used bv lZheng et all (l2007bn . our halo mass function 
produces slightly more halos for M < 2 x lO^^M©, which 
can be corrected by rescaling our halo mass down by 5%, 
and less halos for 2 x lO^^Af© < Af < 3 x W'^^Mq, which 
is possibly due to Poisson noise. 

We derive the corresponding values of M[, A^n-,;,-,, Mq, 
and a in our simulations that closely reproduce the best fit 
nu mber density of cent ral and satellite galaxies provided 
bv lZheng et alj (|2007br ) so that the relative strength of 1- 
halo and 2-halo terms is consistent with the observation 
despite the slight discrepancy in the mass function. In 
detail, we account for the difference in the mass function 
by rescaling M^-^^ and A/q by 5 6%, as our halo 
mass scale in this range is overestimated by 5%, while 
adjusting M[ or a. Note that we do not fit to the observed 
correlation function directly. We show one example of such 
fits (we call 'Model LRG') and the resulting correlation 
function (Tableland Figure [H]). As shown in Table [H 
Af^jjj and Mq for Model LRG are increased by 6% relative 
to the Zheng's values, while cm, M[, and a remain nearly 
the same as Zheng's. 

In Figure [9l the clustering of Model LRG is fairly con- 
sistent with the observed clustering in real and redshift 
space except for r < 0.5h^^ Mpc: the slope of over 
r < 0.5h~^ Mpc is determined by our smallest bin at 
r ~ 0.2h^^ Mpc that is likely subject to the effect of our 
force resolution. In addition, a small discrepancy remains 
near r ~ 2.5h~^ Mpc in real space and r ~ 1 — 2h~^ Mpc 
in redshift space. Again, note that Model LRG is derived 
based on the fitted values of centr al and satellite number 
densities from lZheng et all (j2007b| ) rather than a direct fit 
to the observed clustering data. Although a slight modifi- 
cation in (iVg(M)) of Model LRG may bring the resulting 
^ closer to the observed ^, we find this unnecessary as we 
will focus on quantitative but low-precision comparisons. 

Model LRG produces a shape of ^ which is roughly a 
power law, mainly because of the existence of the shoulder 
in HOD that drives a smaller 1-halo term relative to the 1- 
halo term of dark matter. The shoulder implies a smaller 
fraction of satellites that is mostly from massive halos in 
this case, and so smaller pair counts relative to that o f 
dark matter (see ISeliak|[2000HBerlind fc Weinberg) l2002f) . 
Model 11 in Figure [5] also shows a shoulder in HOD. De- 
spite the small overall satellite fraction, pair counts of this 
model imply a Poisson distribution for AI > IO^^A/q un- 
like Model LRG (§111). Nevertheless, the correlation func- 
tion does not deviate too much from a power-law shape. 
We revisit this and relate this to the effect of evolution in 

m 

We have found that the observed LRG clustering is 
mostly reproducible in our simulations with large values 
of Mi/M rnin and a imposed d irectly at z = 0.3, as moti- 
vated bv IZheng et al.l ( )2007bf ). Meanwhile, the observed 
correlation function of the LRGs does not eliminate the 
possibility of the LRGs having flowed passively from z 1: 
from Figure m those near Model 10 or Model 11 will not 
fare much worse with the observed LRG clustering than 
models near Model LRG would. 



Howe ver, if we take the fit to the LRGs bv lZheng et al.l 
(|20071j ) as concrete, the shoulder in their HOD at z = 0.3 
requires some amount of non-passive fiow evolution, as it 
is difficult to reproduce a large Afi/Afmin and a at z = 0.3 
with passive fiow. In general, we need a process to sup- 
press the LRG satellites in low or intermediate mass halos 
(i.e., 10" Af© - 10^5 Afo) while we need more LRGs in 
very massive halos. As discussed in § 13.11 including dy- 
namical friction will cause our LRG progenitors to merge 
to the center, increasing Afi/Afmin and possibly moving 
a around. This process will be more efficient in lower 
mass halos, as lower mass halos are in general older than 
more massive halos, and so th ey would have had more 
time for this destructiv e process (iKravtsov fc Klvpinlll999l: 



Taffoni et al.ll2003tlZentner fc Bullockll2003l: IZentner et al l 



120051 Ivan den Bosch et al.ll2005bt iTavlor fc Babull l2005f ) . 

For very massive halos, however, we face a contradiction 
of requiring more satellites at z = 0.3 than the conserved 
number density can give, which means we need a new 
source of LRGs in these very massive halos. These new 
LRGs, which are admittedly rare, must have built up enough 
mass between z = I and z — 0.3 through dry merging pref- 
erentially in very massive halos. Note however that some 
models with large initial a (~ 2) do not violate the num- 
ber conservation for high mass halos. For example. Model 
9n4 (§ I3.3[) might reduce to a large Ai^i/Afmin and a after 
non-passive evolution without a necessity for new satellites 
for very massive halos. 

6. GALAXIES ASSIGNED AT Z = 0.3 : FINDING A HIDDEN 
SIGNATURE OF PASSIVE FLOW EVOLUTION 

We have assigned galaxies based on the assertion that 
the number of galaxies is only a function of halo mass 
at the initial redshift. However, the subsequent evolution 
may cause a deviation from this assertion; the cluster- 
ing or the HOD at z = 0.3 may no fonger be a function 
of halo mass alone . Several studies ( jNavar ro et al.l 1997 ■ 
Bullock et al.ll2001l:IWechsler et al.|[20 02: Shc th fc Tormen 
2004HGao et al.l2005tlWechsler et all2005[iZhu et al.l2006 : 



Barker et aLll2006l: ICroton et al.|[2007l : IWetzel et al.l l2007l: 
Jing et al.ll2007t iGao fc Whitd l2007r have indicated that 
details of halo formation history can affect clustering and 
properties of halos of the same final mass. For example, at 
a given final mass, halos that formed earlier show stronger 
concentration, reflecting the high density of the Universe 
at early times, and contain smaller number of satellites 
due to ongoing dynamical friction. It is also shown that, 
for low mass halos, the clustering of early forming halos at 
a gi ven final h alo mas s is strong er than late forming ha- 
los (IGao et aLlboqs.: Wechsler et al.. 2005; .Zhu et al 2006 : 
Barker et al.ll2006l: ICroton et all l2007 r i Wetzel et al.l i2007: 
Jing et al.ll2007l : lGao fc Whitdl2007t) . The clustering trend 
weakens and reverses for very massive halos, in which most 
of the galaxies studied in this paper reside. 

In this section, we look for signatures of passive fiow 
evolution in correlation function at z = 0.3 that cannot be 
parameterized by halo mass at z = 0.3 alone. As a caveat, 
as we assigned galaxies to halos only based on halo mass 
at the initial redshift of z = 1 or 2, we ignored the effect 
of halo assembly history on galaxy occupation prior to 
the initial redshift. It is important to note that the con- 
dition of passive fiow may draw evolutionary signatures 
that are different from non-passive populations. We con- 
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Fig. 10. — The comparisons between § at 2 = 0.3 for passively flowing galaxies (Models 1—31) and ^jj for prompt populations (Models 
IH— 31H) that are assigned at 2 = 0.3 using (Afg(M)} and the second moment of Models 1—31. Left : the ratio of correlation functions in real 
space. Right : redshift space. Solid, dashed, and dotted lines are for three different models, for example, solid line in the top panels are the 
ratio of correlation functions of Model IH to Model 1, long-dashed line : Model 3H to Model 3, and short-dashed line ; Model 5H to Model 
5. The clustering of the two cases is indistinguishable on large scales while different on small scales. 



struct 'prompt' populations at z = 0.3 that have the same 
{N^{M)) and {Ng{Ng — 1))m as those from passively flow- 
ing galaxies but that are randomly distributed among ha- 
los of a given mass. As a result, this 'prompt' population 
is a galaxy population assigned only as a function of halo 
mass at z = 0.3, i.e., without any evolutionary effect. We 
compare the clustering of passively flowing galaxies, i.e., 
evolved populations, with that of the 'prompt' populations 
at z = 0.3. 

In detail, we calculate (A'g(M)) by locating passively 
flowing galaxies (Models 1-11 and Models 21-31) in the 
halos at z = 0.3 and consider the first galaxy in each halo 
as the central galaxy and the rest as satellites. We then 
reassign galaxies back to random halos at z = 0.3 using 
the derived (iVg(M)) for central and satellite galaxies. In 
§ m we claimed that these passively flowing galaxies at 
z = 0.3 are consistent with central galaxies in a near- 
est integer distribution and satellite galaxies in a Poisson 
distribution. We therefore assume these statistics to con- 
struct the 'prompt' populations. The central galaxies are 
located at the most bound particles of halos, while the 
satellite galaxies trace mass. We label these models as 
Models IH-llH and Models 21H~31H, where Models IH- 
IIH are from Models 1-11, and Models 21II-31H are from 
Models 21-31. Some models result in a number of galaxies 
without any associated host halos (about 0.06 ^ 0.2%) at 
z = 0.3, especially from Models 1-11. We simply ignore 
these galaxies. 

Figure [TO] shows the difference in clustering between the 
evolved populations (Models 1-11 and Models 21-31) and 
the prompt populations a.t z = 0.3 (Models IH-llH and 
Models 21H-31H). Both populations produce very similar 
large-scale clustering for the given (A^g(M)), which means 
that it is a good approximation in this scale range to take 
the halo mass as the only variable to decide the number of 
galaxies per halo, at least for populations in very massive 
halos, such as LRGs. The difference remains only on small 
scales, and is larger for models constructed at z = 2 than 



at z = 1. In detail, the clustering of Models 21H-31H is 
larger than Models 21-31 over r = 1 — 2h^^ Mpc, even by 
up to a factor of two for the models with a large number 
of initial satellites (Models 21-23). 

In Figure [51 we have seen the consistency between the 
small-scale clustering and the corresponding pair counts 
among different models. Similarly we will look for any 
remaining difference in the second moments between the 
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evolved populations and the prompt populations, that is, 
the effect of any deviation from a Poisson distribution for 
the evolved satellite populations. Figure [TT] illustrates the 
ratios of (7Vg(M))2 and {Ng{Ng - 1))m between Models 
21H, 25H, 27H and Models 21, 25, and 27, compared with 
the ratio of ^. The ratio of {Ng{Ng — 1))m shows that 
the evolved populations at z = 0.3 arc statistically very 
similar to the prompt populations. The noisy fluctuation 
at the most massive halos contributes only little to the 
overall clustering. As a caveat, pair counts in intermedi- 
ate mass halos with a small satellite fraction are not very 



sensitive to the statistics because they are dominated by 
central galaxy-satellite galaxy pairs, which in turn depends 
on (A^cen(M)) and (A^sat(M)). Because we set these quan- 
tities to be almost identical, the shape of {Ng{Ng — 1))m 
is naturally very similar between the passive and prompt 
populations as long as satellite populations asymptotically 
converge to Poisson statistics. Based on the similarity of 
the resulting pair counts, therefore, the small-scale differ- 
ence observed in ^ appears rooted in the evolutionary effect 
of the halo profiles rather than pair counts. 

From Figure [TUl the difference in ^ is not constant over 
scale: the ratio decreases below 1 for r < lh~^ Mpc but 
increases above 1 near r ^ lh~^ Mpc. The shape of 
the non-monotonic trend implies a broadened ^ for the 
prompt populations (except for Models 9, 10, and 11), 
more close pairs in the evolved populations, and there- 
fore the evolved galaxies taking a steeper radial gradi- 
ent than the prompt populations whose satellite galax- 
ies follow dark matter gradient at z = 0.3. For compar- 
isons, the cluster galaxy radial gradient suggested from 
non-passive studies as well as observations is close to or 



slightly s hallower than mass profile (e.g., Die mand et al. 

"2004: Gao et al.l2004aHvan den Bosch et al.ll2005a:,Nagai fc Kravts 



20051: [Weinberg et al.H20'06r and observational references 



ther ein) with varia tions in morphological type of galax- 
ies (jSpringel et al.l 2001). We try below to provide some 
physical intuition for why the galaxy gradient would be 
steeper than the mass gradient for the passively flowing 
populations. 

Our LRG progenitors in passive flow resided in mas- 
sive halos that correspond to high and rare density peaks 
at early times. At low rcdshift, material from these high 
peaks tends to be more centrally concentrated within a 
host halo than that fro m lower density peaks and the over- 
all m ass distribution (jPiemand et al.l [20051 iMc ' 



et al.l 



I2OO6D . Mass coming from these lower density regions re- 
mains with few galaxies, as we do not allow any new 
sources of galaxies. As a result, a host halo at low red- 
shift will have a high galaxy density in its inner region 
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and little galaxy content accreted into the outer region of 
the halo. In other words, the memory of the typical host- 
halo size at high redshift remains in the process of passive 
flow, resulting in the small-scale difference in ^. On the 
other hand, the trend is reversed for Models 9-11 in Figure 
[TUl where the radial gradient of galaxies appears weaker 
for evolved population, that is, galaxies are on average less 
concentrated than mass in this case. These models started 
with no satellites and probably did not have enough time 
to populate satellites at a close distance while Models 29, 
30, and 31 were able to. In summary, passive flow will 
result in the galaxies being more concentrated than dark 
matter within a host halo, unless these galaxies are recent 
descendants of central galaxies. 

We found that for some cases the difference in small- 
scale clustering due to passive flow evolution is enough to 
cause the correlation functions of evolved galaxies to ap- 
pear more like a power law compared to the prompt popu- 
lations, and vice versa. For example, Models 711, lOH, and 
IIH will give much worse fits to the observed LRG points 
than Models 7, 10, and 11 did in Figure O Therefore the 
effect of evolution could alter the best fit HOD parameters 
to the observed clustering. Note that Model LRG in § [S] 
corresponds to a prompt population. 

The stronger radial gradient for evolved populations could 
also be interpreted in the context of the concentrat ion evo- 
lution of the halos themselves with redshift (iBuUock et al.l 



,— — ,, t pr 

l200lHWechsler et al.l200llZhao et al.l2003allbl ). rather than 
in the context of a radial stratification within a halo with 
redshift. In other words, what we observe could be due to 
the passively flowing galaxies at low redshift and at a given 
final halo mass being more likely distributed among older 
halos, where mass concentration is on average stronger 
(jBullock et al.ll200UrV\^hsler et al.ll2002[ ). as they are ha- 
los from massive progenitors at z = 2 or z = 1. As the host 
halos are fairly massive, it is likely that the corresponding 
difference in large-scale clusterin g is li t tle, as in Figure [TOl 
2005; Wechsler et al.ll2005t iHarker et al.ll2006r 



(iGao et al 



iWetzel et al.ll2G07t IJhIget al.ll2007l : iGao fc Whitell2007l) 

In redshift space, the difference generally appears as a 
relative suppression in the prompt populations, except for 
Models 9, 10, and 11, especially for r < lh~^ Mpc. The 
finger-of-God effect evacuates more small-distance pairs in 
the prompt populations than in evolved populations. As 
the pairwise velocity dispersion within virialized halos in- 
creases with distance (jSheth et al.ll2001[ ). the evolved pop- 
ulations with closer pairs will have a weak e r finge r-of-God 
effect. Also, according to iDiemand et al"] (|2005| ). the ve- 
locity dispersion of material from the rarer peaks is lower 
than that of matter at a given radius from halos, which 
may contribute to the weaker finger-of-God effect as well. 

We consider the possibility that central galaxies in the 
passively flowing populations have not settled down at the 
center of their host halos at z = 0.3 after a series of halo 
merging events. This will also change the small-scale clus- 
tering relative to Models IH-llH and Models 21H-31H 
but probably in an opposite direction to what we have 
observed. We test an extreme case in which no halos 
in prompt populations host central galaxies. In detail, 
we substitute central galaxies in Model 27H with random 
satellite galaxies. Figure [T2] shows that removing the cen- 
tral galaxy strengthens the difference between Model 27II 



and Model 27. In redshift space, the effect is more in- 
tense: the finger-of-God effect is strongly enhanced up to 
r ^ 10h~^ Mpc due to the missing pairs at small separa- 
tion. This implies that the evolved populations at z = 0.3 
probably have a relatively well-positioned central galaxy 
and, again, satellites closer to the center. 

To summarize, the large-scale clustering of populations 
that have evolved through passive flow is reproducible with 
(A'g(M)), and therefore we do not observe significant en- 
vironmental effects on large-scale clustering in this sam- 
ple. The satellite galaxies of passively flowing populations 
closely follow a Poisson distribution, and so they arc indis- 
tinguishable in pair counts from the prompt populations 
with the same (A'g(M)). The effect of passive flow evolu- 
tion, however, appears in the spatial distribution of galax- 
ies within a halo, in that the passively flowing populations, 
unless they are recent descendants of central galaxies, show 
on average more centrally concentrated distribution than 
the prompt populations (i.e., mass profile). 

7. EVOLUTION OF BIAS 

In this section, we study the evolution of bias of pas- 
sively flowing galaxies and compare it to the linear theory 
for passive flow evolution. 

Figure [T3] shows the typical bias evolution for Models 
1-11 and Models 21-31. As the exclusion effect in a halo 
finder is prevailing at initial redshifts, the evolution of bias 
is considered only between z — 0.6 and z — 0.3 for Models 
1-11 (left panels) and between z = 1 and z = 0.3 for Mod- 
els 21-31 (middle panels). The bias factor at z = 0-3 is 
compared to the expected linear bias evolution (jFrvlll996l : 
iTegmark fc Peeble£j|l998l ) where we define 'linear bias evo- 
lution' by unconventionally adopting a scale-dependent 
bias 6(r) and a scale-dependent growth factor of dark mat- 
ter G{r) at z = 0.6 (left, for Models 1-11) or z = 1 (mid- 
dle, for Models 21-31). That is, we define the growth fac- 
tor G{r) at given redshift z as (Cm,z(?')/Cm.zo(''))^^^ where 
zq ^ 0.6 for Models 1-11 and I'for Models 21-31, and 
the bias factor bzg{r) as (Cg(r)/^m(''))^^^ at zq. Then the 
'linear bias evolution' at z expected from zq is derived from 



6iin(r) = (fe.„(r)-l)/G(r) + l. 



(3) 



We then compare the expected bias, &iin(''), with the ac- 
tual bias, b{r) = (Cg('^)/^m(^))^^^ measured at z. 

The bias values at z = 0.3 are consistent with the ex- 
pected bias evolution over r > 2h~^ Mpc in most of the 
models, with only small discrepancies on smaller scales. In 
detail, the bias at z = 0.3 implies that the growth in clus- 
tering becomes slightly impeded over r = 2 — 6h~^ Mpc 
but enhanced over r < 2h^^ Mpc either relative to the 
growth of mass or relative to the linear bias evolution. 
The discrepancies appear bigger for Models 21-31, proba- 
bly due to the longer evolution time and the larger initial 
bias. 

For Models 7n4-10n4 (right panel), the decrease in bias 
is less than that of Models 7-10 due to the smaller initial 
bias values and is also consistent with the expected linear 
bias evolution over r > 2h~^ Mpc. Due to the dominance 
of lower mass halos in these models, their bias factor at 
low redshift deviates from the scale-independent bias on 
smaller scales than Models 7-10. 
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Fig. 14. — The evolution of real-space correlation function of passively flowing galaxies for Models 21, 23, 25, 27, 30, and 31. Correlation 
functions are divided by the correlation function at z = 0.3 (solid lines). Black lines : z = 1. Red : z = 0.6. Blue : z = 0.3. Long-dashed 
lines : the expected correlation functions at z = 0.6 (red) and at z = 0.3 (blue) assuming the linear bias evolution from z = 1 (as defined in 
^ [71. Short-dashed : the expected correlation function at z = 0.3 assuming the linear bias evolution from z = 0.6. One finds that growth of 
clustering is relatively suppressed near r = 2 — 5h~^ Mpc, compared to the growth on other scales. 



To summarize, the evolution of bias for the passively 
flowing galaxies is relatively 'linear' for r > 2h~^ Mpc al- 
though slightly scale-dependent. 

8. EVOLUTION OF CORRELATION FUNCTION: A 
SIGNATURE OF INFALL? 

In this section, we attempt to characterize the evolution 
of clustering of passively flowing galaxies on quasilinear 
scales. 

In Figure [Til real-space correlation functions at z = 1, 
0.6, and 0.3 for a number of models among Models 21-31 
arc divided by a real-space correlation function at z — 0.3. 
The apparent features are the suppression of growth near 
r = 2 — 5h~^ Mpc, and the large growth inside the radius 
(also sec the left panels of Figurc[l]and|4]). The scale where 
the suppression appears is slightly larger than the transi- 
tion from a 1-halo to a 2-halo term where the inflection 
of a power-law model of biased correlation functions oc- 
curs. We compare the growth between z = 1 and z = 0.6 
with what is expected from the linear bias evolution (as 
defined in §[7]) from z = 1: the growth between z = 1 and 
z = 0.3 and between z = 0.6 and z = 0.3 is compared with 
the linear bias evolution from z = 1 and z = 0.6, respec- 
tively. The figure shows that the suppression in growth 
is not explained by linear evolution of clustering of biased 
tracers. 

This effect may be interpreted as pairs at r = 2 — 
5h~^ Mpc moving rapidly to smaller scales as we go to 
lower redshift, that is, an evacuation of pairs on this scale 
as structure grows linearly on larger scales but nonlinearly 
on smaller scales. The feature is less obvious and appears 
at a larger radius for Models 1-11. 

In Figure 1151 we test whether we see the correspond- 
ing feature in the different populations or in the matter 
correlation function. The evolution of matter correlation 
function docs not show an obvious suppression, probably 
because matter is in steady state in flow from a linear to 



nonlinear region of radius. On the other hand, 1-halo and 
2-halo terms of the galaxies do not equally weight the cor- 
responding terms for matter (jSchulz fc Whitd l2006l ) but 
trace the flow of matter differently in the two regions, 
breaking the steady state. 

To examine whether the suppression is a prevalent be- 
havior of biased tracers, we test Models 7n4-10n4. Figure 
1151 shows the suppression in growth of correlation func- 
tion is not obvious at all for Models 7n4-10n4 and so 
the growth of clustering is scale-independent for r > 2 — 
3h~^ Mpc. As a minor point, when compared to the mat- 
ter correlation functions (right panel) , the nonlinear growth 
of the correlation functions for Models 7n4~10n4 happens 
on smaller scales than that of matter; this slightly overpre- 
dicts the linear bias evolution at r = 2 — 4/i^^ Mpc beyond 
a scale- independent growth (in the left panels). The miss- 
ing feature of suppression for Models 7n4-10n4 therefore 
implies that the feature appears preferentially in strongly 
biased tracers (& > 2 at z = 0.3). 

9. CONCLUSIONS 

We have used dissipationless N-body simulations to study 
the effect of passive flow evolution on galaxy clustering and 
halo occupation distributions. We assumed populations of 
progenitor galaxies at z = 1 and z = 2 with a wide range 
of initial HODs, and then studied their properties as they 
flowed passively to z = 0.3. We investigated for a region of 
the parameter space at low redshift constrained by passive 
flow, especially in the halo occupation distribution and in 
galaxy clustering. Our results are summarized as follows. 

Passive flow results in an asymptotic convergence in halo 
occupation distributions and galaxy clustering. The dis- 
tribution of the average number of galaxies per given halo 
mass converges toward a power law for a broad range of 
the initial halo occupation distributions. The values of 
i\/i/A/niin decreases with time, except for the cases with 
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Fig. 15. — The evolution of correlation functions of passively flowing galaxies for Models 7n4 and 10n4(5i,, left panel) and matter (^m> 
right). Correlation functions are divided by the correlation function at z = 0.3 (solid lines). Black lines : z = I. Red : z = 0.6. Blue : z = 0.3. 
Long-dashed lines : the expected correlation functions at z = 0.6 (red) and z = 0.3 (blue) assuming the linear bias evolution from z = 1. 
Short-dashed : the expected correlation functions at z = 0.3 assuming the linear bias evolution from z = 0.6. Cyan lines are at constant 
values and are drawn to help to estimate a scale-dependence of the growth in correlation functions. One finds that the suppression of growth 
near r = 2 — 5/i~^ Mpc is not apparent at all for dark matter or galaxies with small bias {b < 2). 



very small initial Mi/Mmin- The value of a asymptotically 
converges toward unity. Both evolutionary behaviors re- 
sult in (A"g(M)) without a shoulder and thus a shape close 
to a power law at low redshift. 

A similar convergence is observed in the evolution of 
correlation function. While it is not surprising to find 
that the large-scale convergence is consistent w ith the ex- 
pected clustering from h near bias evolution (|Frvl Il996t 
iTegmark fc Peebles! 1 1998h . the intermediate-scale cluster- 
ing also shows a fair degree of convergence in its shape and 
amplitude when the galaxies have evolved from high red- 
shift (z ~ 2). The remaining differences in small-scale clus- 
tering appears consistent with the differences in (iVg(M)). 
Modeling these intermediate and small scales is left for 
future research. 

The asymptotic convergence is enhanced if the galax- 
ies passively flow from higher redshift and so have more 
time to flow. We also find that the satellite populations 
are fairly large at low redshift as a result of passive flow 
evolution. 

We have tried the same tests with galaxy populations 
of a larger number density and found qualitatively similar 
results. The asymptotic convergence in the HODs and 
correlation functions is not an extreme behavior for the 
very high mass halos but is a general result of passive flow 
evolution. We conclude that the convergence of Mi/Mmin 
and a is the result of galaxies asymptotically becoming 
less biased with respect to mass with time. 

Passive flow evolution drives satellite galaxies to con- 
verge toward the Poisson distribution. Other studies for 
non-passive evolution reproduce the Poisson distribution 
as well, implying that the merging, formation, and destruc- 
tion of galaxies happen so that the resulting populations 
remain in the Poisson distribution. 

The comparison of these characteristics of passive flow 
evolution with the observed LRG data hints at any non- 
passive flow processes during the evolution of the LRGs. A 



suggested existence of a shoulder between Mi and Mmin or 
a large a in the HOD of LRGs when fltted to the observed 
clustering dat a (IZheng et al1l2007b . also see Blake et all 
(|2007f) . but see lKulkarni et all (|200<t ) and lHo et all (|2007f) 
for different results) implies that LRGs have not undergone 
a strict passive flow evolution. These discrepancies could 
be due to dry galaxy merging between two LRGs or new 
LRGs arising between the initial and the final redshifts, 
with different efficiencies in different environments. The 
discrepancies should show self-consistent consequences in 
other observational proper ties such as in the luminosity 
funct ion of the LRGs (e. g., Wake et aI1l2006t iBrown et al.l 
I2OOII) . See IWhite et all (|2007l) for an example of such a 
study. 

We compared two populations with the identical (A'^g(M)) 
at z = 0.3, one that is reached after passive flow evolution 
and the other constructed from random halos at z = 0.3, 
to study whether the evolution imposes any distinct sig- 
nature in the second moments of the HOD and cluster- 
ing. The evolved population and the prompt population 
showed no significant difference in the second moments of 
the HOD when both satellite populations are assumed to 
be Poisson at the initial redshift. The effect of evolution 
did not impose a noticeable environmental dependence in 
large-scale clustering. In small-scale clustering, inside ha- 
los, we find that the spatial distribution of the passively 
flowing galaxies is on average more centrally concentrated 
than mass. 

The evolution of bias for passively flowing galaxies were 
consistent with linear bias evolution on quasilinear as well 
as large scales once the scale-dependent growth factor of 
dark matter is considered. 

We identified a suppression of growth of galaxy cluster- 
ing near r = 2 — 5h~^ Mpc. The feature is observed in a 
population with a large initial bias, but not in a population 
with a small bias nor in clustering of matter. 

An interesting extension of this research will be to in- 
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vestigatc the characteristics of passive flow evolution in 
velocity space. Once we understand both spatial and ve- 
locity signatures of passively flowing galaxies, we will have 
a better handle to model the clustering evolution of such 
galaxy populations. We will also extend our work to non- 
passive ly flowing galaxies, using phenomenological galaxy 
merging prescriptions. 



We thank Zheng Zheng for providing his best-fit HOD 
parameter estimates for the LRGs and his valuable com- 
ments on this paper. HS, IZ, and DJE were supported by 
grant AST-0407200 from the National Science Foundation. 
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Table 1 



The initial and final HOD parameters of the galaxies passively flowing from z = 1 



iVioaei 


01 




A/f ^ 
Ml 




f 

,/ ccn 




A'f I A/f 


A/f ^ 


A/f ^ 


f 

7 ccn 








Z = \ 










z = 0.3 






1 


0.5 


2.061 


1.0 


0.4853 


0.58 


0.829 


2.149 


1.847 


0.8592 


0.53 


2 


1 


2.167 


1.0 


0.4614 


0.63 


1.06 


2.414 


1.964 


0.8133 


0.56 


3 


2 


1.929 


1.0 


0.5185 


0.53 


1.67 


2.159 


1.994 


0.9236 


0.48 


4 


0.5 


9.814 


4.0 


0.4076 


0.74 


0.916 


4.223 


3.00 


0.7097 


0.67 


5 


1 


11.16 


4.0 


0.3585 


0.88 


1.11 


5.979 


3.704 


0.6195 


0.79 


6 


2 


11.90 


4.0 


0.3360 


0.97 


1.23 


7.408 


4.281 


0.5779 


0.86 


7 


0.5 


26.50 


10 


0.3774 


0.83 


0.927 


5.390 


3.525 


0.6540 


0.74 


8 


1 


29.43 


10 


0.3398 


0.95 


1.033 


7.353 


4.302 


0.5851 


0.84 


9 


2 


30.39 


10 


0.3291 


0.99 


1.10 


8.646 


4.890 


0.5655 


0.88 


10 








0.3275 


1 


1.04 


9.052 


5.091 


0.5625 


0.88 


11 








0.2141 


1 


1.43 


12.01 


10.39 


0.8651 


0.91 



^Mass of halos are in the unit of lO^M©. 



Note. — We show the input HOD parameters ai z — 1 and the best fit HOD parameters at z — 0.3 
for Models 1-11. Note that Mi and a at z = 0.3 are fitted over M > Mi in order to better describe 
the shape of satellite HODs at the massive end. The value /ccn is the fraction of central galaxies to 
the total number of galaxies. 



Table 2 

The initial and final HOD parameters for galaxies passively flowing from z = 2 



Model 


Q 


Ml/Mmin 


Ml'' 


Mnin'' 


,/ccii 


Q 


A/i/M,„i„ 


Ml'' 




./"ecu 








z = 2 










z = 0.3 






21 


0.5 


1.986 


0.425 


0.2140 


0.60 


0.927 


1.754 


1.601 


0.9131 


0.50 


22 


1 


2.120 


0.425 


0.2005 


0.67 


1.11 


2.226 


1.877 


0.8431 


0.55 


23 


2 


2.144 


0.425 


0.1982 


0.68 


1.46 


2.454 


2.040 


0.8311 


0.56 


24 


0.5 


9.167 


1.7 


0.1855 


0.76 


1.11 


3.162 


2.422 


0.7660 


0.62 


25 


1 


10.23 


1.7 


0.1661 


0.90 


1.15 


4.151 


2.783 


0.6705 


0.72 


26 


2 


10.78 


1.7 


0.1577 


0.98 


1.20 


4.814 


3.022 


0.6277 


0.78 


27 


0.5 


24.34 


4.25 


0.1746 


0.83 


1.131 


3.762 


2.679 


0.7121 


0.67 


28 


1 


26.60 


4.25 


0.1598 


0.96 


1.17 


4.651 


2.978 


0.6403 


0.76 


29 


2 


27.30 


4.25 


0.1557 


0.995 


1.18 


5.134 


3.173 


0.6180 


0.79 


30 








0.1553 


1 


1.18 


5.234 


3.225 


0.6162 


0.79 


31 








0.1056 


1 


1.12 


4.679 


3.735 


0.7983 


0.80 



"Mass of halos are in the unit of IO^^'Mq. 



Note. — We show the input HOD parameters at 2 = 2 and the best fit HOD parameters at z = 0.3 
for Models 21-31. Note that Mi and q at 2: = 0.3 are fitted over M > Mi in order to better describe 
the shape of satellite HODs at the massive end. The value /ccn is the fraction of central galaxies to 
the total number of galaxies. 
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Table 3 

hod parameters at z = 1 with four times the number density of the lrgs 



Model 


a 


Ml/Mmin 


Ml"" 




/ccn 


a 


Ml/Mmin 


Ml'' 




/cen 








z = l 










z = 0.3 






7n4 


0.5 


19.70 


2.6 


0.1320 


0.78 


0.927 


5.502 


1.175 


0.2135 


0.656 


8n4 


1 


22.61 


2.6 


0.1150 


0.907 


0.995 


7.574 


1.403 


0.1852 


0.83 


9n4 


2 


23.91 


2.6 


0.1088 


0.97 


1.132 


9.609 


1.681 


0.1749 


0.83 


10n4 








0.1062 


1.0 


0.952 


10.27 


1.7532 


0.1707 


0.83 



^Mass of halos are in the unit of IO^Mq. 



Note. — We show the input HOD parameters at z = 2 and the best fit HOD parameters at z = 0.3 
for Models 7n4-10n4. Note that Mi and a a.t z — 0.3 are fitted over M > Mi in order to better describe 
the shape of satellite HODs at the massive end. The value /ccn is the fraction of central galaxies to 
the total number of galaxies. 



Table 4 

The five-parameter HOD at z = 0.3 for Model LRG 



Model 


a 


M(=' 


M' ■ " 


/ccn 


LRG 


1.86 (1.84 


6.875 (6.875) 


0.8226 (0.7764) 


0.931 (0.932) 



"Mass of halos are in the unit of 10^" M©. 

'"We use Mo = 3.209 (3.056) x lO'-'M© and ctm = 0.556 (0.556). 

"Num bers in the parentheses are the best fit values by IZheng et al.l 
(|2007bl ) for our fiducial cosmology. 

'^The number density of galaxies for Model LRG is 9.817 x 
W~^h^ Mpc""^. This is only slightly different from the fiducial number 
density 10~*/i"^ Mpc~^ that we adopted for passive flow. 



Note. — HOD param eters at z = . 3 in ou r simulations that correspond 
to the best fit HOD bv lZheng et all (|2007bD for the observed LRG clus- 
tering. For comparison, when Model LRG is fitted to the three-parameter 
HOD (eq. [T]), we find a = 1.86, Mi/M^in = 12.2, Mi = 6.790 x 10"Mq, 
and Mmin = 0.5563 x lO^Mg. 



